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1. Introduction 

The interplay between disorder and many-body interactions continues to be one of 
the main topics of condensed matter physics. In one- dimensional quantum many-body 
systems without disorder, the Bethe ansatz (BA) has proven to be a valuable tool, giving 
access to the energy spectra and eigenstates of certain so-called integrable models |^. 
At first glance, it appears that the integrability, namely the existence of infinitely many 
integrals of motion, seems to preclude any applications of the method to disordered 
systems. However, in 1984 it was shown ^ how to apply the BA method to the 
Kondo problem |^ of a single magnetic impurity in a bath of conduction electrons. 
Further developments led to the construction and solution of integrable spin chains 
with embedded spin defects [|], 0. 

A different approach to integrable impurity models was considered in references [§] 
where impurity vertices are introduced by varying the local interaction parameters while 



preserving integrability. These studies have stimulated further investigations |T0[ of 
such impurities in various systems. The resulting models have impurity terms with 
couple to the charge degrees of freedom, and look fairly similar to generic impurity 
terms. However, the energy spectrum is independent of the spatial distribution of the 
defects, and there is no localization of the ground-state wave function, unlike what is 
expected for generic impurities. This peculiar behaviour can be understood by the fact 



that integrability implies a purely forward-scattering mechanism at the impurities ITp 
There is no reflection and thus no possibility of destructive quantum interference which 
could lead to a localization. 

Back-scattering can be introduced into integrable models by choosing suitable 
boundary conditions (BC). Sklyanin |jTlI] proposed a systematic approach to construct 



and solve integrable quantum spin systems with open BC. Central to his method are the 
so-called reflection equations (RE) |12] which are the boundary analogues of the Yang- 



Baxter equations (YBE) [T^. Together, the YBE and RE imply the integrability of a 
model which can then be constructed as usual by the algebraic approach of the quantum 
inverse scattering method (QISM) [|14|. The finite-size corrections of the corresponding 



energy spectra and the asymptotic behaviour of correlation functions follow predictions 
based on boundary conformal field theory [|1^]. We remark that the BC of Sklyanin 



are called "open" in order to distinguish them from the more often used periodic and 
the free BC. Although the term "open" seems to suggest particular transmission and 
reflection properties, this is not necessarily implied. The combination of open BC and 
integrable impurities has been considered in references [ll6|, p!7|, p^, p^. Of particular 



interest is the case where the forward-scattering impurity is directly coupled to a back- 
scattering open boundary |TB| , P^ . This combination may lead to physically relevant. 



yet completely integrable models. 

In the present article, we construct two kinds of integrable impurities for a fermionic 
small-polaron model with general open BC. Due to the fermionic nature of the model, 
we employ the graded version of the QISM ||2^, For well-separated impurity vertices 



Integrable impurities for an open fermion chain 



3 



located within the bulk, the local interaction terms involve the two neighbouring sites as 
usual ip. Placing the forward-scattering impurities at the back-scattering boundaries, 
we derive a Hamiltonian with rather general boundary terms which may be interpreted 
as sources and sinks of particles at the boundaries. Using the graded YBE and the 
graded RE, we derive the BA equations, and obtain expressions for the eigenvalues for 
special cases of the Hamiltonian. In addition, we discuss the ground-state properties in 
the thermodynamic limit. 

The paper is organized as follows. In section ^ we introduce the small-polaron 
model with general open EC. In section ||, a class of integrable impurities is constructed 
by shifting the spectral parameters of local Lax operators at arbitrary sites in the bulk. 
By embedding the impurity fermion vertex at each boundary of the model, we construct 
a class of integrable impurities with perfect back-scattering in section |[ In section H, we 
study the algebraic BA solutions for those impurity models. The ground-state properties 
are discussed in section ||. Section |^ is devoted to a discussion and conclusion. 



2. The small-polaron model 



We consider the small-polaron model [^, which describes the motion of an additional 
electron in a polar crystal. The Hamiltonian reads 

N N N 



+ p+rij^ + a+ajy + + p_n^ + a_a\ + /S-a^, (2.1) 

where J is proportional to the overlap integral, V denotes the electron-phonon coupling 
and W is the energy of the polaron. The boundary coefficients p±, a± and /9± are 
Grassmann variables, with p± even and a±, (3± odd. Hermiticity of the Hamiltonian 
requires = (3± and = P±- The fermionic creation and annihilation operators a] 
and flj satisfy the usual anticommutation relations 



{aj,a^} = {a],al} = 0, {aj,al} 



6. 



(2.2) 



and rij = cLjCLj- The i?-matrix and local monodromy matrix are explicitly given as 



R 



12[Ul,U2) 



/a'l 





\ 


-tb'L c" 







c" ibl 







\ 


-a'L 


) 



[2.3] 



and 



Lj{u) 



b' 



b',)n. 



ca- 



-ic'a] 



(a'_ + ib'_)n- 



[2A) 



respectively. They satisfy the graded Yang-Baxter algebra (YBA) 



Rl2{Ul,U2) T{Ui) T{U2) = T{U2) T{Ui)Ri2{Ul,U2) 



[2.5) 
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where 



T{u) = Ln{u)---L2{u)Li{u), 


(2.6) 


and 

1 2 






(2.7) 


Here, ®5 is the supertensor product 






(2.8) 


with the parity -P(l) = 0, P{2) = 1 such that the i?-matrix correspond to the "null" 


parity case P(a) + P(/?) + P(7) + P(5) = 20 . We parametrize the coupling 


' parameters 


J, V and W as 




J =1, 


(2.9a) 


V = -2c2(0), 


(2.96) 


W = 2s2(0)tan(w) + 2c2(0). 


(2.9c) 


The entries of the P-matrix (|2.3|) and the monodromy matrix (|2.4|) are 




ai =e^'(wl)ef Ms2(mi-m2), 


(2.10a) 


=e?K)e^'MsoK-M2), 


(2.10&) 


C" = c' = S2(0), 


(2.10c) 




{2.10d) 




(2.10e) 


Cq{u±w) 

Co{u)co{w) 


(2.10i) 


where we introduced the convenient notation 




Sn{u) = sin('U + nrj), Cn{u) = cos{u + nrf). 


(2.11) 


Throughout the paper, we therefore use rj and w for the parametrization of the model 


parameters V , and W . 




In a previous article |24| , we proved that the model is integrable under the 


conditions that the boundary K± supermatrices 








-'"'^(3 3) 


(2.12) 


take the form 






(2.13a) 




(2.136) 


K^2 = - ""-"'fi't^ [es2{u) - e^s.2{u)] [eM^) - e^-2(«)] , 

<+'^-S2(0) 


(2.13c) 


K2I = ^^2(0) [^-^2(^i) ^+S-2{U)\ [^+S2{U) ^_S-2{U)\ , 


{2.13d) 




(2.13e) 
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K+ = ^.S2{u + [eM^) - esoiu)] , (2.13/) 

= - ""^'ffi'f^^ [es.iu) - csoiu)] [e^^) - e-soiu)] , (2.13,) 

Kt, = - ^-^'".^tlnr^''^ [^-^^(") - - • (2.13/.) 



Here we would like to emphasize that although we ||2J] construct the general boundary 
ii"-matrices ( p.l2| ) for the small-polaron model (|2.1|) by the Lax pair formulation, we did 
not figure out the form of the RE corresponding to more general boundary i^'-matrices 
( p.l2| ). In the above expressions, we further defined 

p± = S2(0)cotV^±. (2.14) 

The parameters tp± control the strength of the boundary potential terms, whereas a± 
and j3± in ( |2.13aD -( p.l3/j| ) are the parameters characterizing the fermion sources and 



sinks at the boundaries. The Hamiltonian ( |2.1| ) can be obtained as usual as an invariant 



of the commuting family of transfer matrices t{u) 

t{u) = StTo[K+{u)T{u)K4u)T-\-u)] (2.15) 
by taking the derivative at a special value of the spectral parameter u. Namely, 

(2.16) 



2/fT(0) + Stro ( ^K+{u] 

u=0 



du 



u=0 



with Stro denoting the supertrace with respect to the auxiliary space. 
3. Integrable impurities in the bulk 

In this section, we construct integrable impurities which appear in the bulk part for 
the fermionic small-polaron model with general open BC. If the quantum i?-matrix of a 
fermionic system has the difference property of spectral parameters, the associated Lax 
operator with additional parameter also satisfies the graded YBA, i.e., 

12 2 1 

Rl2{Ul-U2) LiUi + u) L{u2 + v) = L{u2 + v) X (Mi + I/)i?i2 (Mi - ^2) . (3. 1) 

Therefore one can construct a class of integrable impurities for the fermion model with 
both open and periodic BC by shifting the spectral parameters of local monodromy 
matrices at arbitrary sites in the bulk. The associated monodromy matrix is given as 

T{u) = Ln{u) ■ ■ ■ L.m{u + u.m) ■ ■ ■ Li{u) (3.2a) 

T-\-u) = L^\-u) ■ ■ ■ L^\-u + vm)--- L-^\-u), (3.2&) 

where the parameter i/^ characterizes the impurity strength at site m. Now we suppose 



that the supermatrices K± are the solutions of the graded RE [|25|, ^ 

1 2 
Ri2{ui - U2) K- {ui)R2i{ui + U2) {U2) = 

2 1 

K- {u2)Ri2{ui + U2) K~ {ui)R2i{ui - U2), (3.3a) 
Rll^'^'^'i-ui + U2) k^^'{u,)Ri2{-ui -U2- 4r/) kf\u2) = 

kf'{u2)R2i{-ui -U2~ Ar]) i:^ti(^^)/?swst.(_^^ ^ ^^^^ (3_3^) 
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where R2i{u) = VRi2{u)V with the graded permutation 

= (3.4) 

and iStc denotes the inverse operation of the supertransposition Sta in the space a. The 
i?-matrix possesses the unitary property Ri2{u)R2i{—u) = p{u)I, with p{u) a scalar 
function of u. It follows that the double-row transfer matrix (|2.15|) may be considered 
as a generating function of the infinite conserved quantities such that the i^'-matrices 
are the solutions to the RE (|3.3o|) and ( p. 3^) . From the relation ( p.l6|) , it is not difficult 
to obtain the Hamiltonian for the open fermion chain with an impurity located at site 



TV 



H 



J2 H,,^, + -L,{0)KUO)LY\0) + 



StTo[K^{0)L'j,{0)L-\0)] 



i=2 

jj^m,m—l 



Stro[i^+(0)] 



m~-2-^m,m— 1 



where 



A 

-'^m,m—l 



x,_i(o)x;.(o)z-i(o)x-_\(o), 

X„_i(0) X„(z/„) 1.-1(0) i„i_i(0), 
X„_i(0) X„(0)i-i(i/™)X„i_i(0), 

x„_i(o) L'^{v^) L-^\Q) 



(3.5) 

(3.6a) 
(3.66) 
(3.6c) 
(3.6c/) 



The prime denotes the derivative with respect to the spectral parameter u. The 
interactions of the open fermion chain with the impurity are shown schematically in 
figure 0. 

In order to simplify the algebraic calculation for the construction of such an 
integrable impurity for the model ( |2.1| ), we let w = 0. The supermatrices ( p . 1 2| ) become 
(up to a normalization) 

1 f—Soiu — ip^) a_So(2u) 
(3_Sq{2u) Sq{u + 
S2{u-'ip+) a+S4(2u) 

/?+S4(2m) S2{u + iIj+) 



K_{u) 
KAu) 



so(^-) 



(3.7a) 
(3.7&) 



Then we check that the supermatrices (|3.7q| ) and ( p.76| ) satisfy the graded RE ( |3.3g| ) 
and ( |3.36| ), respectively. From ( p.5| ), after some algebra, we obtain the Hamiltonian for 
the small-polaron model with both general open BC and an integrable impurity located 
at site m as 



H 



N 

i=2 



(b) 



H 



(b) 
N 



m,m—l,m—2 ' •• 



m,m—l,m—2 



m.m— l.m— 2 



with 



a\a,_^ + a\_^a, + 2c2(0) ri,n,_^ - 2c2(0) n 



(b) 



^2(0) 



Co{iIj-) — ia_ a\ — if3_ 



(3.8) 

(3.9a) 
(3.96) 
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m-1 



-e — e- 



m-2 



N-1 



N 



Figure 1. Graphical representation of the interactions in the chain with the impurity 
located at an arbitrary site m together with boundary impurities. 



^(b) 



^2(0) 



H 



(h) 

■m,m— l,m— 2 



S2(0) - S2(i^m) 



Co(^/'+) - ia+ aj^ - a 



N 



(d) 

m,m— l,m— 2 



a 



^ a 

m—2 m 



+ h.c. 



H. 



(c) 



m,m— l,m— 2 



go(t^m)g4(0) 

A(z/„) 

2C2{0)sl{Um 



[^m^m-2 



^m-2^m-l + 



^m^m-2'^m-l 



''^m-2'^m'^m-l 



A(z/„ 



-n 



m-l'^m-2"'»n 



h.c, 



where 



A(n) 



and h.c. denotes the hermitian conjugate with {ip^ 



h[^^ and H^' are the general boundary terms which are responsible for the sources 
and sinks with particle injection and ejection at the boundaries; -ffm,m-i,m-2 consists 
of nearest- and next-nearest-neighbour hopping terms involving the sites m, m — 1 and 
m — 2; h':^^_^ contains an onsite potential and density- density interaction terms 
between neighbours and next-nearest neighbours; and Hm\m-i,m-2 involves current- 
density interactions (see figure 0). The Hamiltonian in the presence of more than one 
impurity can easily be constructed, if the two nearest impurities are still well separated. 
In this case, the Hamiltonian reduces to a sum over the isolated impurities like in the 
case of the Heisenberg periodic chain [p|, p|, p!0 . 



(3.9c) 



(3.9cO 

(3.9e) 



(3.9/) 



(3.10) 
V. Here, 



r(b) 



4. Integrable impurities coupled to the boundaries 

Kondo-like impurities of local impurity spins coupled to ID strongly correlated 
conduction electrons have attracted much interest [0, especially in the context 
of the BA solution [^. To every complex- valued i^'-matrix solution of the RE ( 3.3c| ) 
and (|3.3fe|) , one may construct a class of "regular" solutions [0, i.e., 

k_{u) = L{u)K_{u)L-\-u), k+{u) = K+{u) (4.1) 
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to the same RE. In order to study Kondo impurities for ID electron systems |T9|, |27 
it is better to construct "non-regular" , i.e. non-trivial operator-valued, solutions, to the 
RE. Wang and coworkers [0 constructed a class of integrable impurities coupled to 
each boundary of the spin-^ Heisenberg XXZ chain by a special choice of boundary K±- 
matrices, i.e., K± = 1. But in their approach the parameters characterizing the strength 
of the magnetic impurities — related to our potential impurities via the customary 
Jordan- Wigner transformation [|l| — disappear in the Hamiltonian as well as in the BA 
equations. Here we present a different approach to integrable impurities: from "regular" 
solutions of the graded RE (|3.3q| ) and (p. 3^ ), we construct a class of integrable impurities 
^ [lO| coupled to each of the boundaries of a fermion chain with general open BC. 
We stress that these impurities are not Kondo-like. If we embed two fermionic impurity 
vertices at the boundaries, 



T{u) = Lr{u + VrjLNyu) ■ ■ ■ 
T-\-u) = Li\-u + ue)L^\-u) 

one can show that 

U4u) = T{u)K4u)T-\-u) 



u) ■ ■ ■ Li(u)Le{u + Ui 



l-^\-u)l;H 



-U + Vr 



(4.2a) 
(4.26) 



(4.3) 

also satisfies ( |3.3o| ) and so does the solution Li{u + Ui)K^(u)Lj^{—u + Vi). It follows 
that there exists a family of transfer matrices 

t{u) = StTo[KMU-{u)] (4.4) 

and its members commute with each other for different spectral parameters. Similarly to 
( p.l6| ), we can formulate the explicit expression of the Hamiltonian for an open fermion 
chain with boundary impurities. 



H 



N 



i=2 



1 



{StTo[K^{0)L'^{iyr)L;\iyr)] 
+StTo[K+{0)Lr{iyr)L'j,{0)Lj,\0)L;\ur)]} 



+ -L,{0)Le{iy,)KL{0)Lf{ue)Li\0) + L^{0)L',{u,)Lj\iy,)Li\0). 



(4.5) 



Substituting ( p.7«|) and ( |3.76| ) into 



the corresponding Hamiltonian is given as 



H 



N 

+ 



i.i-i 



^2(0) 



^2(0) 



A(M,)so(^^ 



(i) 

N,r 



A(z/,)so(^-) 



(4.6) 



where 



H 



(b) 
N 



C2(^^r)S-2(^^r)So(^^ 



n 



Hj:^^ = S2(0) [s2{ip+) n,. - is2{iyr)a+ a]- + is^2{i^r)(^+ O-r] 



(4.7a) 
(4.76) 



H 



(i) 

N,r 



- ^2(0)02(0) [so{il)+) n, - 2iso{ur){a+ - /3+ a,.)] un 
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-e — e — e — e- 



12 N-1 N 

Figure 2. Impurities coupled to each of the boundaries. 



— e — e — e — e — e — ^ 

2 m-2 m-1 m N-1 

Figure 3. Integrable impurities situated at the boundaries. 



- C2(0) S2(0)so(^+) Un + 2isl{vr){a+ a]v + 13+ a^v) 

- S2(0) So(z/r + ^+) a-jv^r - ■So('^r - aJoAT , 



(b). 



H 



(b) 
(i) 



H 



(i) 



(4.7c) 

(4.7d) 
(4.7e) 
(4.7/) 



l,r ^ £,tp+ ^ tp-,a+ P-,K -t^e), 

where i?^*^) can be interpreted as fermion sources and sinks with particle injection 
and ejection at the boundaries and at the impurity sites. However, unhke the previous 
Hamiltonian ( |3.8| ), the boundary parameters and impurity parameters are both involved. 
iJ*-'-' describes the interaction between impurities and boundaries (see figure 0). 

On the other hand, if we move the impurity in the bulk to each boundary of the 
chain as shown in figure ^ with the monodromy matrix 

T{u) = Ln{u + un) ■ ■ ■ Lmiu) ■ ■ ■ Li{u + ui), (4.8a) 
u] = i^i ^(-M + z/i ) ■ ■■l^,^^[-u) ■ --L^^^-u + z/iv), (4.86) 

) apart from the numbering 

r ^ iV, ^ - 1, £ 1, 1 ^ 2. 

Although the eigenvalues of the open chain do not depend on the position of the 
impurities in the bulk due to the absence of back-scattering, the open boundary is 
a perfect back-scatterer with vanishing transmission at each end of the open chain for 
a± = /?± = 0. Moreover, it is easy to obtain a model with the impurities coupled to 
each boundary together with / well separated impurities (see figure H) at positions 
for i = 1, . . . ,/, i.e. 

N 



T-\-u) = L^\-u + I.,) ■ ■ ■ L-^ 
one finds that the Hamiltonian is same as 



^2(0) 



J=2 



+ 



A(z/r)so(^+) 



(i) 
N,r 
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m-5 m-1 

^ — - - — e — - - — 

m-6 m-4 m-3 m-2 m N r 

Figure 4. Two well-separated bulk impurities at sites mi — m and m2 = m — 4 
together with the boundary impurities. 

/ 

' / j \ mi,mi~l,mi~-2 ' mi,mi — l,mi—2 ' mi,mi — l/m^— 2/ ' \ 



The terms H^^\ H^^ are the same as in 

(|]7^-(|TJ), and the terms i?^^) , i?^'^) 



and H^^^ are given in ( p.9tj) -( ^.9j ). To keep the Hamihonians ( p.8|) , (|4.6| ) and 
hermitian, the parameters rj and ip± must be purely imaginary, u real and a^. = P±. 
All terms in the Hamiltonians are needed to ensure the integrability of the models. 
In the next section we shall proceed with the algebraic solutions for the small-polaron 
model with different kinds of impurities in the most interesting special case of perfectly 
back-scattering boundaries without sources and sinks. 

5. Bethe ansatz solution for finite chains 

Following the method of [jrT], we shall study the algebraic BA solutions for the open 
fermion chain with different kinds of impurities. We first note that the general open 
BC spoil the pseudo vacuum state. Therefore it seems difficult to solve the models with 
general open BC by means of the QISM. We thus restrict ourselves to the simpler 
situation a± = 0, (3^ = in the following. In this case, the Hamiltonians (|3.8|) , 
( [4.6|) and ([4.9|) do not contain any Grassmannian source and sink terms and are charge 
conserving. They still contain the potential impurities and are perfect back-scatterers. 
Thus these Hamiltonians are ideal for the proposed investigation of the interplay of 
forward-scattering bulk impurities with backward-scattering boundaries. 

Let us for simplicity first consider the Hamiltonian ([4.6|) . In the case a± = 0, 
(3± = 0, the Hamiltonian ( |4.6|) comprises 

/fj?) = [s2(0)c2(0)so(V'+) - slii^Mi^-,)] n^, (5.1a) 
Hj>^ = [s2(0)c2(0)so(^+) + sl{0)co{^+)] n„ (5.15) 

H'^,r=-S2{0) So(z/r + ^+) O-jv^r - So(l^r - V'-h) O-Jo-Af - S4(0)So(V^+)n,.njv (5.1c) 
rib) rrib) 



and H['\ Hf\ and follow as in ( |4.7c| )-( |4.7J ). H^^"* describes the boundary 



impurities, contains the interaction terms with exchange coupling between the 

bulk and the impurities (see figure As mentioned before, this Hamiltonian conserves 
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the particle number due to the absence of sources and sinks with particle injection and 
ejection at the boundaries. 

Now we proceed to establish the Bethe eigenvectors for the Hamiltonian ([4 .61) with 
( p.laD -( p.lcj) by means of the algebraic BA [Q. Let 



Tiu) 



A B 
C D 



-u 



A B 
C D 



(5.2) 



be the monodromy matrices acting on the pseudo vacuum state defined by cij |0) 
0, j = 1, ■ ■ ■ , A^. Then we have 

A|0) = s^(m)so(u + z/^)so(m + m,)|0), 



D|0) 



U)S2{U + VI,)S2{U + Z/r)|0), 



5|0) =0, 
C|0) ^0, 



A\Q) 



5 [T{-u)] 

(-l)^S^ {U)S2{U - yi)s2{u - Ur) 



6 [T{~u)] 



|0), 



|0), 



5|0) =0, 
C\0) ^ 0, 

where the quantum determinant 
define 

A B 



is 6 [T{-u)] = A^{u)A{u - i^e)A{u - 



Uju) 



C D 



From (|4.3| ), we then have 

1 



A 
D 



1 



[-so{u - ip.)AA + so{u + ip-)BC], 
[-so{u - ^^)CB + so{u + ij^)DD]. 



Noting the following form of the graded YBA 

T-\-u)Rui2u) Tiu) =T{u)Ri2{2u) T'^ 
it is possible to derive the commutation relation 

^2(0) 



BC 



S2{2u) 



(DD-AA). 



With the help of the graded RE ( 3.3q|) we obtain — after a lengthy calculation 
commutation relations 

Aiu)Civ) = '4^^^^P^Civ)Aiu) - ^^(0)^^(2^) ciu)Aiv) 



(5.3a) 
(5.35) 
(5.3c) 
{5.3d) 

(5.3e) 

(5.3/) 

(5.3^) 
(5.3/i) 

Let us 
(5.4) 

(5.5a) 
(5.5&) 

(5.6) 

(5.7) 
the 



So{u - V)S2{U + v) 

S2{0)so{2v)si{2u) 



+ 



S2{2v)s2{2u)s2{u + v) 



S2{2u)so{u — V 
C{u)D{v), 



(5.8a) 
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So{u - V)S2{U + v) 

^2(0) . 

S2{u + v) ^ ^ ^ ^' 
where we introduced the transformation 

S2{2U) 



so{u - v)s2{2v) 



{5.8b) 



(5.9) 



From (|5.3G| )-( |5.3/i| ) and ( [5 .41 ), we can choose an M-particle excitation as 

Mv, ■ ■ ■ vm)) = C{vi) ■ ■ ■ C{vMm. (5.10) 
Using the commutation relations (|5.8q|) and ( ^.8^ ), one obtains the eigenvalue A of the 



(5.11) 



transfer matrix 

r(ti) |$(fi ■ --Vm)) = A{u;vi ■ ■ ■t;A/)|$(fi ■ ■ -fju)), 



where 

A{u;vi ■ --Vm) 



■1) 



AT 



S2{2u)5[T{-u)]so{^,) 
X <( S2{u - ip+)s2{u - iJ-)so{u - ue)so{u + ue)so{2u)sl^{u) 

, , , ^ TT S4(m + Vo)s2(m - t^,) 

]X S^[u-Vj)S2{u^V^) 

+So(m + i\)j^)sq{u + '\\)S)S2{u + V()S2{u - V^Si^{2u)S^^ {u) 

.^SQ{u\Vj)S.2{u-Vj) 

XS2{U + Ur)s2{u - l^r) \ \ 7 ' V 

f_l So{u-Vj)s2{u + Vj) 



(5.12) 



provided that 



„y^S-l(^i + '^m)5-l(wj-Z/J AJ-^ S2{Vj+Vk)S2{yj-Vk) 



M 



S-.2{Vj + Vk)s^2{Vj - Vk) 



(5.13) 



for all j = 1, . . . , M. In the above BA equations, we have shifted the parameter 
Vj Vj — rj. From the relation ( 2.16| ), the eigenvalue E of the Hamiltonian 
for a± = P± = follows as 



E 



MO) 



cot ^_ + cot ip^ + 2{N + 1) cot 2r7 - 2 coi{vr - 2r]) - 2 cot{iye - 2r]) 

M 



2s2(0) 



(5.14) 



On the right hand side of (|5.14| ) we have dropped a multiplicative term l/2cos(2r7) 
as is customary [|ll|]. We emphasize that the spatial position of the impurities in the 
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chain, although clearly important in the construction of the Hamiltonian, is irrelevant 
for the BA equations and the ground state energy 0, ^, This is the mathematical 
formulation of the physical absence of backscattering for these impurities [1^ . 

The BA equations and the eigenvalues for the Hamiltonian ( |3.8| ) can be obtained 
similarly as, 

S-i{Vj + llj-)S-i{Vj + ip+)sf''^~'^\vj) 



for all j 



. . , M, and 



M 

n 

k = l 



S-2{Vj + Vk)S-2{Vj - Vk) 
S2iVj + Vk)s2(Vj - Vk) 



(5.15) 



E 



-S2(0) 



M 



cot + cot ifjj^ + 2N cot 2?7 — 2 cot(z/m — Sr^) — 



2s2(0) 



(5.16) 



For the Hamiltonian 



we find 



siivj - 4)-)si[Vj - i}+)s\ 



n 



SljVj + l^m)Sl{Vj - Vm) 



M 

n 

fc = 1 



S-2{Vj + Vk)S-2{Vj - Vk) 
S2{Vj+Vk)s2{Vj-Vk) 



(5.17) 



where 11™=! denotes the product over the / isolated impurities in the bulk as well as 
the boundary impurities. The energy spectrum is given as 

cot ^„ + cot + 2(A^ - / + 1) cot 2?7 - 2 ^ cot(z/^ - 27]) 



E 



-S2(0) 



m=l 



M 



-E 



2s2(0) 



S^i{Vj)si{Vj) 



(5.18) 



When comparing ( ^7] ), (§]TB|) to ( pl3[ ), ( P^ and ( |^1^ , (pl^ , we see that using 
( p. 17] ), ( p. 181 ) corresponding to the Hamiltonian ([4.9|), we can reproduce the results for 
the other two Hamiltonians (|3.8|) and (|4.6|) . Namely, with / = 1 and z/^ = z/^ = we get 
the result for (pj) with (|5T5|) , (|5l6D and with / = and z/^, z/; ^ we find for (U) 
the BA equations ( ^.131) with energy (|5.14|) . Thus the Hamiltonian ([4.9|) contains the 
other two Hamiltonians as special cases, although the construction by QISM proceeds 
independently. We note that care has to be paid to the varying number of sites N when 
doing this procedure. 



6. Ground-state properties in the thermodynamic limit 



We note that the bulk terms of the Hamiltonians (|3.8| ), (^4.6| ) and ([4.9|) are equivalent 
to the counter part of the ID Heisenberg XXZ model with periodic BC via a Jordan- 
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Wigner transformation. The finite-size corrections and thermodynamics for the XXZ 
model with or without boundary magnetic fields have been studied in many papers 
p9| , pO| . As mentioned before, the Jordan- Wigner transformation does not preserve the 
boundary terms nor the impurity terms due to its nonlocality. The BA equations we 
obtained provide a more meaningful description of the ground-state properties due to 
the presence of the boundary potential terms and the impurity parameters. In what 
follows, we shall study the ground-state properties for the resulting models following 
the scheme in g |18|, |9|, |0|. 

For convenience, let us first redefine the variable Vj ivj. Then, taking the 
logarithm, we rewrite the BA equations ( ^.13| ) for the Hamiltonian ( [4.6|) as follows: 

2TTlj = 2N9{vj, 7]) + 9{vj,'i{j+ - v) + ^("^j) ip- -rj) + 6{vj + z/^, rj) + 9{vj - z/^, r]) 

M 

+9{vj + ue, 7]) + 9{vj - ue, v) - ^ ^(^j - ^fe^ '^v) + (^{^j + Vk, 2r]), (6.1) 

k = l 



for all j = 1, ■ ■ ■ , M, where the two-body phase shift p, [T^, [1^ is 



smhivj + irf) 

9[Vj, ri) = im — = 2arccot (tanhf ,• cot r/j . (6.2) 

sinh(t>j — ir]) 

We now define V-k as —Vk and define vq = 0. Then the density of the roots {vj} can be 
defined as 

dZN{v) Ij 
PNiv) = ^^, Z, = - (6.3) 



and the finite-size BA equation (|6.1| ) becomes 

1 1 ^ I 

(^{^^v) + i^[e^Hv) + 9^'^\v)\-— e{v-VkM\. (6.4) 

k=-M 




where 

9^'\v) = 9{v + Ur,ri) + 9{v - Ur,ri) + 9{v + ue,ri) + 9{v - ue,ri), (6.5a) 
^(b)(^,) = 0(^y^ 27]) + 9{2v, 2ri) + 9{v, tp+ - r]) + 9{v, i). - rj). {6.5b) 

We note that the first two terms in ( |6.5^) arise due to the non-periodicity of the chain, 
whereas the last two terms are due to the boundary potentials at sites 1 and N. Taking 
the thermodynamic limit and differentiating (|6.4| ) with respect to the spectral parameter 
V, we have 

Pooiv) = - h'iv, v) + 7^ k'^'(^) + ^^'^^'(^)l 1 - r dupMd'iv - 2r7), (6.6) 



n \ ' ' " 2N 



27T 



-A 



where the integration boundary A is determined by 



Poo{v)dv = ^^ + 0{N-'). (6.7) 

■A 
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The prime denotes the derivative with respect to v. Due to the hnearity of ( |6.6| ), one 
may formally solve the following three linear integral equations: 

P^^\v) = -0'iv,v)-^ r dup^^\u)9'iv-u,2v), (6.8a) 
p«(^) = _ 1 /■ dup^^{u)e\v - 2r/), (6.86) 

p(^)(t;) = V)'(t;)-^ !\upt\u)e\v-u,2r^), (6.8c) 

In this way, the solution of ( |6.6|) can be expressed as 

Poo{v) = p^^{v) + ^[pt{v) + pt\v)i (6.9) 

where p^^{v), j^p'i^{v) and 2^p£^(f) are the contributions of the bulk, the impurities 
and the boundary effect to the root density, respectively. The ground-state energy 
( ^.14| ) is minimized at the cutoff A in the thermodynamic limit as discussed in [p9| , |3C 



Following the argument in we find the cutoff A = oo such that the particle 

density is M/N = 1/2. 

By using the Fourier transforms, the formal solutions to the equations (|6.8Q|) - (|6.8(j ) 

read 

p^K'ri (6.10) 

where 

_ /"OO 

e{uj,T]) = / e'{v,ri)e''^''dv. (6.11) 



From the residue theorem, we obtain 



pSH^) = (6.12a) 

7] cosh 

.^-^ 4cosh77-f COshTT^Z/m 

p('\v) = V ^ , (6.126) 

°° ^ ?7 cosh ^{v + Vrn) cosh ^{v- I/m) ' 



m=r,< 



1 / 

Pt\v) = - [pt^\u) + p'^^^u)] e—du, (6.12c) 



(6.13a) 
(6.136) 



where 

.(be)/ N _ 2 sinh(f - 2r])uj + 4 cos sinh(f - t])uj 

^'^^ ~ sinh f + sinh(f - 2r])uj 

-fhnW \ 2sinh(f + 7] - ?/^+)a; + 2sinh(f + ?7 - V^.)!:^ 

^ ' sinhfcj + sinh(f -2r/)cj 

Here p^^\uj) and p"^\uj) are the contributions to the root density from the boundary 
effect and the boundary potential terms, respectively, due to ( |6.56| ). Then, from ( |5.14| ), 
we also obtain the ground-state energy in the thermodynamic limit as 

4 sin^ 2n 
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with 



Eq 



-P+ - P- 



2{N + l)cos2r] + 2sm2ri cot{um - 2r]) . (6.15) 



The boundary energy is given by 

2 sin^ 2r] 



E. 



dv- 



cosh 2v — cosh 2ri 
p+ — p_ — 2 cos 2?7 + 2 sin 2ri cot( 



(6.16) 



2r^). 



We thus note that the boundary potential terms do not only enter the expression for 
the ground-state energy explicitly as — p+ —p-, but also implicitly via of ( |6.12t| ). 



We remark that in 



boundary magnetic field terms, corresponding to ( p.l3fc| ), 
do not contribute to the root density due to the lack of free boundary parameters in 
the boundary ii'-t-matrices. The presence of boundary potentials (magnetic fields) and 
the impurity parameters changes the asymptotic behaviour of the BA equations ( |6.1| ) 
resulting in string solutions different from those discussed in the papers [^. Indeed, 
either the boundary parameters p± (or or the impurity strength parameters Um, 
Ui and Ur affect the boundary string solutions to the BA equations. It is obvious that 



the ground-state energy of the bulk is same as in the periodic case p9|. In general, the 
boundary states are associated with complex roots of the BA equations. 

Analogously, we obtain the ground-state energy (|6.14| ) for the Hamiltonians (|3.8| ) 
and ( [4.9|) . The differences to the ground states for these Hamiltonians are only the 
contributions from the impurities expressed in p^{v). Thus we get for the Hamiltonian 



4 cosh cosh 1^1/^ 



T] cosh ^{v + u^) cosh ^{v-u^) r] cosh ' 
— p+ — p_ — 2A^ cos 2ri + 2 sin2r7 cot(z/m — 2?]). 
For the most general Hamiltonian ([4.9|) , we have 



P^(v) = E 



4cosh|^f cosh|^i/„ 



Eq 



T] cosh ^{v + Um) cosh ^{v -iy„ 

- p+ - P- - 2{N + 1 - f) cos 2r] 

f,r,e 

+ 2 sin 2?7 cot(z/, 



/ 



rj cosh 



2r/). 



(6.17) 
(6.18) 

(6.19) 
(6.20) 



m=l 



Further thermodynamic properties such as compressibilities and susceptibilities can be 
calculated as demonstrated previously in M, |10|, . Results will be presented elsewhere. 



7. Conclusions and discussion 



In the present work, we have considered the interplay of integrable impurities and general 
open boundary conditions for the example of the small-polaron model. The impurities 
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have been constructed via inhomogeneous shifts in the spectral parameters of the Lax 
operators such that the YBE are satisfied. The boundary terms are taken to obey 
the RE. In both cases, we dealt with the graded version of the equations due to the 
fermionic nature of the particles and the boundary terms. Thus by construction, the 
model remains integrable. We have shown that this is true even when placing the 
impurities directly at the boundaries. 

The most general boundary terms considered in (|2.1| ) include fermionic particle 
sources and sink terms as well as more standard density terms. However, these 
linear terms in creation and annihilation operators contain coefficients that are 
odd Grassmann variables. Thus a straightforward physical interpretation appears 
problematic. Representing these coefficients a-t, (3± as additional fermionic operators 
a-t, aj_, we arrive at a chain with two additional sites but without sources and sinks. 

The boundary terms coupling to the particle density can be viewed as potential 



impurities — much like in the Anderson model of localization |3T[] — situated at the 
boundaries. For the special case with only these potential impurities and the integrable 
impurities present, we solve the BA equations and compute the ground-state energy in 
the thermodynamic limit for half-filling. We find that the solution is possible regardless 
whether the integrable impurities are located within the bulk or at the boundaries. 

The two types of impurities enter the expressions for the ground-state energy 
additively. Thus the simultaneous presence of both purely forward scattering integrable 
impurities and purely refiecting boundary potential terms does not seem to change the 
physics in a substantial way. We therefore do not expect to see the onset of localization 
as might be expected from the form of the boundary impurities. 

An analogue of the integrable impurities can be found in the case of light waves. 
Consider a long strip of glass, interspersed with pieces of bifringent material of the same 
index of refraction as the glass. Then as circular polarized light enters the strip and 
reaches the first bifringent slab, its plane of polarization will be rotated by an angle 
6 oc all where /i represents the length of the first slab and a is the material specific 
rotation angle per unit length [^. There is no refiection at the contact due to the 
identical indices of refraction. After the next slab, we have 6 = ali + al2 and so on. 
Thus the net effect of the bifringent slabs is a rotation of the plane of polarization, i.e., 
a change in the overall phase of the wave function of light just as for the electronic wave 
function in Ref. P, |^, |10|- In this picture the integrable boundaries satisfying the RE 
then simply correspond to perfect mirrors at both ends of the strip. 
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